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Based on a unified theoretical treatment of the ID Bogoliubov-de Genes equations, the super- 
fluidity phenomenon of the Bose-Einstein condensates (BEG) loaded into trapped optical lattice 
is studied. Within the perturbation regime, an all-analytical frameworli is presented enabling a 
straightforward phenomenological mapping of the collective excitation and oscillation character of 
a trapped BEG where the available experimental configurations also fit. 
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INTRODUCTION 

Harmonically trapped Bose-Einstein condensates 
(BECs) offer a great chance to understand the macro- 
scopic quantum phenomena such as phase coherence [l|-|3| 
and matter wave diffraction [J] . Condensates loaded in a 
periodic potential forming an optical lattice (OL) [5|, |6| 
may show a rich dynamic picture of oscillations as Bloch 
oscillations [7|, |8[ , Belieav and Landau damping [91, llO| , 
Landau- Zeeman tunneling [7[, and the a ppea rance of the 
superfluid oscillation of condensates [ll|, \\^ . Although 
schemes of controlling the dynamics of BECs have been 
profusely described [8|, Il3l4l5 1 , the underlying physics of 
some of these studies appears hidden under numerical 
analysis. The description in terms of excited states or 
Goldstone modes not only highlights the main cause of 
this behavior, but also enables the characterization of 
the BEC dynamics in universal terms. Thus, we select 
a platform of the BEC loaded simultaneously into a 
harmonic traps and an optical lattice to characterize 
the phenomenon of superfluidity as well as the dynam- 
ical properties and tackled the problem analytically. 
The purpose of this letter is to derive a perturbative 
treatment which allows explicit closed solutions for the 
phonon dispersion relation and to reveal the effect of 
BEC configurations on the superfluidity phenomenon. 
The method has undergone the test of comparison with 
experimental evidences with success. 

Systems such as the cigar-shaped trap schematically 
represented in the upper panel of Fig. [1] can be con- 
sidered quasi-one-dimensional (ID) confinement. Within 
the framework of mean field theory, the physical char- 



acteristics of a BEC loaded in such trapping profile are 
ruled by the time dependent nonlinear Gross-Pitaevskii 
equation (GPE) 0. 
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where \id is the self-interaction parameter, m is the al 

/27r 
kaline atom mass, Vext{x) = ^muj^x^ — Vj, cos^ I — 

represents the harmonic trap potential plus the periodic 
potential caused by the counter-propagating lasers with 
Vl the laser intensity, d its laser wavelength, and wq the 
frequency of the harmonic trap. For Eq. ([1]), it is pos- 
sible to prove rigorously , the existence of ground states 
for any Aid and Vl. Moreover, the set of ground states is 
orbitally stable and the ground states have a Gaussian- 
like exponential asymptotic behavior for any Aid, regard- 
less of Vl value [l7[ ■ An important consequence of this 
mathematical fact is the stability of those physical mag- 
nitudes, such as no explosion and no damping as function 
of time, which are described by operators defined in the 
Hilbert space of the ID GPE dJ). This result is par- 
ticularly related to the superfluidity properties, among 
others physical phenomena, of the harmonically confined 
condensates loaded in optical lattices. 

The collective excitations, or so-called Goldstone 
modes of the BEC, can be obtained by applying a small 
deviation from the stationary solutions |^o) of Eq. ([1]), 



|*(t)) = exp(-i/it/fi) [|*o) + |u)exp(-iwt) 
-I- \u) exp(— iwt) + \v*) exp(iwt)] , 
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FIG. 1: (Color online), (a) The oscillating BEC loaded into a 
ID optical lattice within a parabolic trap (grid) . The elliptical 
spot symbolizes the condensate with a center of mass, X{t), 
oscillating around the trap bottom, (b) Collective excitation 
energies hujk of the first 10 modes calculated as function of 
dimensionless parameters A = Xid / (lofvjJo) for Vo = , and 
(c) as function of Vo = VL/huoo for A = 2 and d/lo = 0.25. 
The solid (dashed) lines represent odd (even) modes. 



which corresponds to Unearizing the time-dependent non- 
Unear Schodinger equation in terms of amplitudes \u) and 
\v*), fi being the chemical potential and w the mode or 
phonon frequencies [18|. Inserting |^(t)) into Eq. ^, we 
obtain the Bogoliubov-de Genncs equations 
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where £ =p2/2m + Vea,t-^ + 2Ai£) (*o |*o)- Three cou- 
pled nonlinear equations for |^o)j \u) and \v) must be 
solved simultaneously, which incorporate the harmonic 



trap and the stationary optical potential. In general, 
this is a very onerous task and it is not always possible to 
extract transparent solutions giving reliable information 
on the BEC dynamics. If the harmonic trap potential is 
switched off from Eqs. dD) and ^ it is possible analyti- 
cally extract reliable information of the ID condensate in 
a periodic potential as the Bloch oscillation and stability 
of the solution [19|, |20[. This corresponds to the ho- 
mogeneous case where the phonon wavevector q = q^ex 
is a good quantum number and the excited frequency 
Wq = w((7) is a continuous function of q. For wq 7^ 
the wavevector q is no longer a good quantum number 
(inhomogeneous case) and the system (IT]-[2]) provides a 
set of discrete excited states ujk labeled by k = 1,2, .... 
By assuming a weakly-interacting Bose gas and not too 
strong laser intensities, the self-induced nonlinear inter- 
action and the OL potential can be considered as pertur- 
bations with respect to the harmonic trap potential. Ac- 
cordingly, compact solutions for ^ and the spatial shape 
of the order parameter |\I'o) , as determine d by relevant 
parameters of the condensate are obtained [21 1 . 



NORMAL MODES 



Considering the nonlinear term Xio (^0 l^o) a-nd the 
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periodical potential Vl cos 



as small terms com- 



pared to the confined harmonic trap potential strength, 
the solutions of the system ^ can be cast in terms 
of the complete set of harmonic oscillator wave func- 
tions {\ipn)} ^- \u) = J2n=o ^n \ipn) and \v) = 
^j^_Q Bn IV'n) and the coefficients An and i?„ arc ex- 
panded in form of series An = An + An + ...; Bn = 
Si^^ + B^'' + ..., where the quantities aI'^ and B^P - A' 
, V^. Using Eqs. ^ we get 



E^i"(^i") 
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-Vo (V'fcl cos2aa; |i/'„) + 2A {^Jk |*o> (*o| i'n) 
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with a = 2Trlo/d, Iq = y^h/mujo, A = AiD/(Zofia;o), Vo = 
Vi/fiajn. Taking advantage of the procedure developed in 
Ref . |2l| and solving simultaneously the system (HJ , it is 
possible to show that the independent phonon frequencies 
LOk are given by 



Wo 



fcH 

Vo 



A 

/2^ 



-1 



2r(fc + i/2) 



0FA: 



2 cxp(-a2)[i,(2«2)_l] 



/2^ 



exp I —a 



2 2 

^ 2 ^ 2 



4(a) 
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where Lk{z), r(z), Ei(z), Chi(z), and C are the Laguerre 
polynomials, the gamma function, the exponential inte- 
gral, the cosine hyperbolic integral, the Euler's constant, 
respectively. Finally, jk are numbers. Sk and pk being ex- 
plicit functions the dimensionlcss parameter a |23l . [24 1 . 

In Fig. [T] the analytical solutions for the frequencies 
LOk arc graphically represented for the first 10 modes. 
Panel (b) displays Wfc for attractive (A < 0) and repulsive 
(A > 0) cases at Vq = 0. Also in panel (c) the influence 
of the periodical potential on the phonon modes being 
checked for fixed values d/lo = 0.25 and A = 2. In the 
first case, the collective oscillations show an almost flat 
dispersion as a function A. Note that the mode for A: = 1 
has the frequency value wq of the harmonic trap |25| . In 
panel (c) is also seen a blue-shift renormalization of Wfc 
can be noted due to the presence of the OL. 

The normalized eigenvectors |$fe) = (I'^fc) : l^fc)) can 
be cast as 



^ H'k 
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with/fc^„ = (-1)*'' 

^_^^(fc-m)/2^,^2a) 

expf— a^), ip(.z) being the Laguerre polynomials, p = 
{k+m—\k — m\)/2 and m+k = even number. Notice that 
the quasiparticle amplitudes jufc) and |v^) present parity 
inversion symmetry property, i.e. if the index k is an even 
or an odd number we are in presence of two independent 
subspaces where the wave functions |$fc) become sym- 
metric (even mode) or antisymmetric (odd mode) with 
respect the transformation x — > —x.. 

According to Eqs. (0) and (jH]), for a given time f, 
the probability density |(a;| ^fe(t))| of the excited states 
along the x -axis shows oscillations with well defined max- 
ima which are quenched according to the exponential be- 
havior exp (— x^/Zq) . The position of the maxima and 
minima of the axial density |(x| ^fe(i))| are linked to 
the minima or maxima of the combined potential func- 
tion VjhLtjQ = 0.5(a;/Zo)2 — T^ cos^ (2-^x11^. Finally, an- 
other important limit is reached when the optical lattice 
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FIG. 2: (Color online) (a) Universal spectral power contribu- 
tions to the BEC oscillations within the perturbation regime. 
The inset shows the mode number, fcm, where the maximum 
value of the spectral power is attained, (b) Center-of-mass 
position. Experimental data from Ref. |Tl| are represented 
by open squares and full circles, while by solid lines show the 
calculations after Eqs. ((S)) and ([7}. 



is turned off. In this case the standard solutions (in a 
perturbative sense) of Bogoliubov-de Gennes equations 
for the inhomogeneous case are directly obtained from 
Eqs. © and © taking Vo = . 



SUPERFLUIDITY. 

To characterize the dynamics of a BEC loaded 
in an OL we evaluated the expectation value of 
the center-of-mass position X(t) —l^(t)\x\^(t)), with 

oo 

1^(0) == E ^fe l<^fc(i)) • Here, the set of eigensolutions 

fc=0 

is chosen as |(/>o(i)) = exp(— i/ii/fi) |^o) , I0fc(i)) = 
exp(— z/zi/fi) [exp(— iwfet) |Mfe) + exp(ia;fci) \v*^\ for fc 7^ 0, 
and the coefficients {Cfe} are obtained under certain ini- 
tial condition. In our case, we consider that, at i = 0, 
the OL is absent and condensate is located out of equi- 
librium at certain distance Xq from the origin, i.e. the 
order parameter |^o) is centered at a; = Xq with an ex- 



pectation value X{0) = Xq. At, t > the system may or 
may not be loaded into the OL periodical potential. 

A straightforward calculation, by keeping terms up to 
first-order in A and Vq , yields that the dynamics of the 
center-of-mass is ruled by the equation X{t) = Xo(i) + 
XA,Vo(.'t): where 



Xoit) = Xoexp^ Fk cos {ujk+i-uJk)t, (7) 

fc=0 

Fk = cxp{-X^/2ll)-X'^''/[{2ll)'' kl] is the spectral power 
and X\y„ (t) is a linear function of A andVo, with negligi- 
ble contribution to X(t) for typical experimental setups. 
In Eq. [7]the condition ujk=o = is used. 

The spectral power contribution to these oscillations, 
in terms of the mode number, fc, exclusively dependent 
on the relative initial position, Xq/Iq, as displayed in 
Fig. [2] (a). For a given initial displacement, the mode 
k = km with maximal contribution to the Xo{t) is 
given by the equation In (^X^/2ll') — Hk^ — C = 0, with 
Hn = X]fe=i ^/^- Thus, one can see the larger the rel- 
ative initial displacement, Xq/Iq, the higher the main 
contributing modes, as shown in the inset of Fig. [2] (a), 
and wider their diffusion. This has energetic implications 
since the mode frequencies, uj^, obtained from Eq. ^ are 
tunable with the nonlinear interaction strength and the 
OL parameters. Without OL the vibrational level spac- 
ing Awfc = LOk+i — ujk ~ wn which make the expectation 
value Xo(i) = Xq cos wpi |25| . 

Figure [21(b) shows the superfiuid oscillation for a BEG 
of ^^ Rb extracted from Ref. 11 measured in a static 



Rb extracted from Ref. 
magnetic trap with and without a ID periodic potential, 
as simulated by our approach. Using Eqs. ([5]) and (O we 
are able to reproduce the reported experimental center- 
of-mass position as function of time without using any 
fitting parameter. The oscillations observed in Fig. [2] (b) 
correspond to vibrational level spacing Aoj^ of the two 
independent subspaces with even and odd modes that 
make X{t) oscillate with a frequency near the harmonic 
value ujf). The small frequency shift of the condensate 
loaded into OL, and observed in Fig. [5] (b), is directly 
linked to the renormalization of the atomic mass of the 
system moving in a periodical potential. 

In summary, we presented a unified analytical descrip- 
tion for the collective excitations (phonon frequencies ujk, 
Eq. dS]), the excited states wavefunctions |$fc) , Eq. ©), 
the ID superfluidity oscillation and the dynamics (center- 



of-mass position Xo(i) Eq. ([7])) of BEG systems loaded 
in an optical lattice. 
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where ^ is restricted to m-\- k = even number, p = (/c + 



